home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / claqhp.z / claqhp
Text File  |  1996-03-14  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAQQQQHHHHPPPP((((3333FFFF))))                                                          CCCCLLLLAAAAQQQQHHHHPPPP((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAQHP - equilibrate a Hermitian matrix A using the scaling factors in
  10.      the vector S
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CLAQHP( UPLO, N, AP, S, SCOND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED, UPLO
  16.  
  17.          INTEGER        N
  18.  
  19.          REAL           AMAX, SCOND
  20.  
  21.          REAL           S( * )
  22.  
  23.          COMPLEX        AP( * )
  24.  
  25. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  26.      CLAQHP equilibrates a Hermitian matrix A using the scaling factors in the
  27.      vector S.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      UPLO    (input) CHARACTER*1
  32.              Specifies whether the upper or lower triangular part of the
  33.              Hermitian matrix A is stored.  = 'U':  Upper triangular
  34.              = 'L':  Lower triangular
  35.  
  36.      N       (input) INTEGER
  37.              The order of the matrix A.  N >= 0.
  38.  
  39.      AP      (input/output) COMPLEX array, dimension (N*(N+1)/2)
  40.              On entry, the upper or lower triangle of the Hermitian matrix A,
  41.              packed columnwise in a linear array.  The j-th column of A is
  42.              stored in the array AP as follows:  if UPLO = 'U', AP(i + (j-
  43.              1)*j/2) = A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-
  44.              j)/2) = A(i,j) for j<=i<=n.
  45.  
  46.              On exit, the equilibrated matrix:  diag(S) * A * diag(S), in the
  47.              same storage format as A.
  48.  
  49.      S       (input) REAL array, dimension (N)
  50.              The scale factors for A.
  51.  
  52.      SCOND   (input) REAL
  53.              Ratio of the smallest S(i) to the largest S(i).
  54.  
  55.      AMAX    (input) REAL
  56.              Absolute value of largest matrix entry.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAQQQQHHHHPPPP((((3333FFFF))))                                                          CCCCLLLLAAAAQQQQHHHHPPPP((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      EQUED   (output) CHARACTER*1
  75.              Specifies whether or not equilibration was done.  = 'N':  No
  76.              equilibration.
  77.              = 'Y':  Equilibration was done, i.e., A has been replaced by
  78.              diag(S) * A * diag(S).
  79.  
  80. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  81.      THRESH is a threshold value used to decide if scaling should be done
  82.      based on the ratio of the scaling factors.  If SCOND < THRESH, scaling is
  83.      done.
  84.  
  85.      LARGE and SMALL are threshold values used to decide if scaling should be
  86.      done based on the absolute size of the largest matrix element.  If AMAX >
  87.      LARGE or AMAX < SMALL, scaling is done.
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.